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Abstract 

The Higgs triplet model (HTM) predicts the existence of a pair of doubly 
charged Higgs bosons i^ ±± . Single production of i/ ±=t via cy collision at the next 
generation e + e~ International Linear Collider (ILC) and the Large Hadron elec- 
tron Collider (LHeC) is considered. The numerical results show that the production 
cross sections are very sensitive to the neutrino oscillation parameters. Their values 
for the inverted hierarchy mass spectrum are larger than those for the normal hier- 
archy mass spectrum at these two kinds of collider experiments. With reasonable 
values of the relevant free parameters, the possible signals of the doubly charged 
Higgs bosons predicted by the HTM might be detected in future ILC experiments. 
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1. Introduction 

During the past decade, neutrino oscillation experiments have provided us with very 
convincing evidence that neutrinos are massive particles mixing with each other [1]. This 
exciting breakthrough gives a strong motivation for physics beyond the standard model 
(SM), since the SM itself only contains three massless neutrinos [2]. Various ways to 
go beyond the SM have been proposed in order to accomodate this observation. The 
Higgs triplet model (HTM) [3, 4] is a simple extension of the SM with an SU(2) L - 
triplet Higgs boson of hypercharge Y = 2 whose vacuum expectation value (VEV) va 
provides Majorana neutrino masses without introducing right-handed neutrinos. 

The HTM [3, 4] generates a Majorana neutrino mass via the product of a triplet 
Yukawa coupling hij and a triplet VEV v^. This model predicts the existence of seven 
physical Higgs bosons which are two CP-even neutral bosons h° and H°, a CP-odd 
neutral one A , a pair of single charged bosons H ± , and a pair of doubly charged Higgs 
bosons H ±:iz . These new particles can produce rich phenomenology at present and in fu- 
ture high energy collider experiments. Since charge conservation prevents doubly charged 
scalars from decaying to a pair of quarks, iJ ±± carry the lepton number violation and 
would give rise to a distinctive experimental signal, which has lower background. This 
fact has lead to many studies involving the doubly charged Higgs bosons H ±±L in the lit- 
erature. For example, it has been shown that such particles can be produced with sizable 
rates at hadron colliders via the processes qq — > H ++ H [5, 6] and qq' — > _P ±± _P T [5, 7, 
8]. Using these production mechanisms, discovery of if ±± at the Large Hadron Collider 
(LHC) has been studied in Refs. [9, 10]. 

There is a direct connection between the Yukawa coupling and the neutrino mass 
matrix which gives rise to phenomenological prediction for processes which depend on 
[11]. The branching ratios of the lepton flavor violating (LFV) decay processes U — > Ijhh 
and li — > Ij'-f contributed by the doubly charged Higgs bosons if ±± are strong depend on 

in the context of the HTM. Recently, it has been shown [12] that the current upper 
experimental limits for some of these LFV decay processes can give serve constraints 
on the parameter space of in the HTM. Taking into account these constraints, in 
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this paper, we consider single production of if ±=t via e^ collision at the next generation 
e + e~ International Linear Collider (ILC) [13, 14] and the Large Hadron electron Collider 
(LHeC) [15]. Although single production of the doubly charged Higgs bosons via cy 
collision has been studied in Refs. [16], in this paper, we will focus our attention on 
the effects of the neutrino oscillation parameters on its production cross section and see 
whether the neutrino mass hierarchy can be disentangled via this kind of production 
processes at the ILC and LHeC. Our numerical results show that the production cross 
sections at these two kinds of high energy colliders are very sensitive to the neutrino 
oscillation parameters. The production cross section of the process cy — > e + H for the 
inverted hierarchy mass spectrum is significantly larger than that for the normal hierarchy 
mass spectrum. In the sizable range of the parameter space of the HTM, their values 
can be large enough to be detected in the future ILC experiments. 

Our paper is organized as follows. The essential features of the HTM are briefly re- 
viewed in section 2. In this section, we also summarized the constraints on the HTM from 
some LFV decay processes. The effective production cross sections for the subprocess 
ej — > l + H at the ILC and LHeC are calculated in sections 3 and 4. Our conclusions 
and discussions are given in section 5. 
2. The essential features of the HTM 

The Higgs triplet model (HTM) [3, 4] is one of the appealing scenarios for generating 
neutrino mosses without the introduction of a right-handed neutrino. In this model, a 
complex SU(2)l triplet scalar with the hypercharge Y = 2 is introduced to the SM, 
which can be parameterized by 



A 



1 A+/V2 A++ ^ 



v A° -A+/V2 j 
A nonzero VEV value < A > gives rise to the following mass matrix for neutrinos 



(1) 



rrii 



Mj = 2h tj < A >= V2hijU A . (2) 
The symmetric complex matrix hij = e,fj,,r) is the Yukawa coupling strength. The 



interaction of the triplet scalar with lepton doublet Lj = (v iL , l iL ) T is given by 



£ 



-h l3 [vJ L Cu jL A» - -j=(i&Cl jL + lJ L Cu jL )A + - l* L Cl jL A ++ ] + h.c, (3) 



here C is the charge conjugation operator. 

The neutrino mass matrix m^ can be diagonalized by the Maki-Nakagawa-Sakata 
(MNS) matrix Vmns [17]. Then the couplings can be written as following 
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Where mi, m.2, HI3 are the absolute masses of the three neutrinos. The phase 5 is the 
Dirac CP-violating phase, ipi and </? 2 are referred to as the Majorana phases [4, 18] and 
there are < ipi, ip 2 < 2tt. s^- = sin^-, c^- = cos^- with < ^ < 7r/2, in which 0y are 
the mixing angles. The explicit expression forms for have been given in Refs. [8, 9]. 

After the eletroweak symmetry breaking, there are seven physical massive Higgs 
bosons (if ±=t , H^ 1 , h° , H° , A ) left in the spectrum for the HTM. The doubly charged 
Higgs bosons if ±=t are entirely composed of the triplet scalars A ±:t , while the remaining 
eigen states are mixtures of the doublet and triplet scalars. Such mixing is proportional 
to the triplet VEV which is generally small even if z/ A assumes its largest value of a 
few GeV. Thus, H ± , H°, A are predominantly composed of the triplet scalars, while the 
SM-like Higgs boson h° is mainly composed of the doublet scalars. 

It is well known that the doubly charged Higgs bosons H ±i: have contributions to the 
LFV decays k — > ^7, lj — > Ijhk [19]. In the HTM, their branching ratios depend on the 
neutrino mass matrix parameters, the VEV value ua, and the mass m H ±± of the doubly 
charged Higgs bosons [12]. The current experimental upper limits for some of these LFV 



decay processes can give severe constraints on the relevant free parameters. In the HTM , 
the LFV decay process /i — > cy proceed via exchange of iJ ±± and H^ 1 at one-loop. In 
the case of ignoring the electron mass in the final state and all lepton masses in the loop, 
the branching ratio Br(/i — > cy) can be written as [12] 

Briery)* , (6) 

where a = 1/137 is the fine structure constant, G F = 1.17 x !0~ 5 GeV~ 2 is the Fermi 

constant. In above equation, it has been assumed m H ±± = m H ± = m#. The current 
experimental upper limit for the LFV process /j, — > is Br exp (fi — > ey) < 1.2 x 1CT 11 
[20], then we have 

H 2.6 x 10-(^)'. (7) 

Using Eq.(4), there is 

2u 2 A (hh + ) ij = m\5ij + [V MNS diag(0, Am^, ^rnl 1 )V^ NS ] ij , (8) 

where Am^ = m\ — m\ and Am^ = m\ — m\. Since the sign of Am^ is undetermined 
at present, distinct patterns for the neutrino mass hierarchy are possible. Am^ > is 
referred to as normal hierarchy (NH) with m\ < m2 < m% and Am^ < is defined as 
inverted hierarchy (IH) with m 3 < m x < m 2 . Combing Eq.(7) and Eq.(8), there is the 
following relation 

<lv\ > 1.96 x 10 A ( 2 °° GeV ) 2 [V M Nsdtag(0, Am 2 21 , Am^y^]^. (9) 
"'// 

In our following numerical estimation, we will consider this constraint on the triplet 
VEV v A to calculate the cross sections for single production of the doubly charged Biggs 
boson H via cy collision at the ILC and LBeC. 

3. Single production of the doubly charged Higgs boson H at the ILC 

The LHC can generate very massive new particles and will essentially enlarge the 
possibilities of testing for new physics effects, while the ILC is also need to complement 
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the probe of the new particles with detailed measurement. For the ILC, the center of 
mass (cm.) energy y/s = 0.5 — ITeV and the typical integrated luminosity £i n t = 0.5 — 1 
ab^ 1 is currently being designed [13, 14]. An unique feature of the ILC is that it can 
be transformed to 77 or e7 collision with the photon beam generated by laser-scatting 
method. The effective luminosity and energy of the 77 and ey collisions are expected to 
be comparable to those of the ILC. In some scenarios, they are the best instrument for 
discovery of the new physics signatures. In this section, we will consider single production 
of the doubly charged Higgs boson H predicted by the HTM via e7 collision. 
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Figure 1: The Feynman diagrams for the partonic process e {jp%) + 7(^2) ~~ > 
l + (p 3 ) + H—(p 4 ) in the HTM. 



From discussions given in section 2, we can see that the doubly charged Higgs boson 
H can be produced via e7 collision associated with a positive lepton. The relevant 
Feynman diagrams are shown in Fig.l, in which / + denotes the lepton e + ,/x + or r + . For 
the process e~~(pi) + 7(^2) — > ^ + (P3) + H (^4), the renormalization amplitude can be 
written as 

M = M a + M b + M c = — 2eKl - v T (p 3 )C- l P L ( J / 1 + ^)Yu( Pl )eM 

(Pi +P2) 2 



4:eh e i 
(j>2 ~ Pa) 2 - m 2 H 

2eh e i „.T/\n/i 1 



v T (jpz)C- l P L u{p x )(2 Pi - p 2 )%(p 2 ) 



(/;| ni) y (Psh^i- ^)C- L P L u( Pl )e,(p 2 ). (10) 

Here Pl = (1— 7s)/2 is the right-handed projection operator and £^(^2) is the polarization 
vector of the photon 7. The explicit expressions for the Yukawa coupling constants h ee , 
h efl , and h eT can be written as [8, 9] 

hee = ^^[m 1 e-^44 + m 2 ^ 2 Ci3 + ^ (25 -^ ) )^3, (H) 
V2z/ A 



[mie llpl c 12 c 13 (-s 12 c 23 - e lS c 12 s 13 s 23 ) + m 2 s 12 c 13 (c 12 c 23 



e~ tS s 12 s 13 s 23 ) + m^'^ si 3 c l3 s 23 



;i2) 



y/2; 



[mie lLpl c 12 Ci 3 (s l2 s 23 - e l5 c 12 c 23 s 13 ) + m 2 c 13 s 12 (-c 12 s 23 - 



e- l5 s 12 s 13 c 23 ) + m 3 e l ^-^s 13 c 13 c 23 ] 



(13) 
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Figure 2: At the ILC, the production cross section <j\ for the subprocess ey — > e + H 

as a function m for different values of the m# in the NH(a) and IH(b) cases. 



After calculating the cross section a(s) for the subprocess — > l + H , the effective 
cross section <J\(s) at the ILC with the cm. energy -y/i = ITeV can be obtained by 
folding the cross section d(s) with the photon distribution function / 7 / e (a;) [21] 

r0.83 

= / . , S(s)f y/e (x)dx, (14) 

JDlj/S 

where x — s/s, in which ^ is the cm. energy of the subprocess e~7 — > l + H . In above 
equation, we have neglected the mass of the lepton l + due to m# >> m ; +. 

It is obvious that the production cross section o\ is dependent on the neutrino os- 
cillation parameters, the triplet VEV h>&, and the mass parameter m#. Direct searches 
for H ±:L have been performed at LEP [22], Tevatron [23], and HERA [24] via the pro- 
duction mechanisms e + e~ — > H ++ H and qq — > H ++ H . Both of these production 
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mechanisms are depend on only one unknown parameter m#. The lower mass limits for 
ma from the Tevatron searches are mg > HOGeV ~ 150GeV [23]. In this paper, we take 
rriH as a free parameter and assume that its value is in the range of 300GeV ~ 600Ge^. 
The current experimental upper limit for the LFV process /i — > has give severe con- 
straint on the triplet VEV as shown in Eq.(9). In our numerical estimation, we take 
the minimal value given by Eq.(9). In this case, the production cross section o\ is only 
dependent on m# and the neutrino oscillation parameters. 

There are nine neutrino oscillation parameters: three neutrino masses m^ (or m^), 
three mixing angles 9ij, the Dirac phase S and two Majorana phases tpx and (fi2- According 
the current constraints on the neutrinos and mixing parameters from neutrino oscillation 
experiments [25], we use the following values in this article 

Am 2 21 = m\-m\~ 7.59 ± 0.20 x 10 _5 e1/ 2 , 

y -0.69 J 

Am^ = m 2 3 -ml = -2.36 ±0.11(±0.37) x 10~ 3 eV 2 for IH, 

Am 2 31 = m\ -m\ = 2.46 ± 0.12(±0.37) x 10~ 3 eV 2 for NH, (15) 




9 12 = 34.4 ±1.0 , 023 = 42.8^ 



+10.7 
-7.3 



13 = 5.6+|°(< 12.5)°. 



(16) 



Information on the mass mo of the lightest neutrino, the Dirac and Majorana phases can 
not be obtained from neutrino oscillation experiments. To simply our paper, we will take 
5 = and (pi = <pi = to calculate the cross sections for single production of H in 
the cases of NH (mo = mi < m<i < m^) and IH (mo = < m\ < 777.2). For the free 
parameter mo, we will assume its value in the range of O.OOleV < mo < leV. 
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Figure 4: Same as Fig. 2 but for the subprocess ey — > t + H 



From above discussions we can see that the production cross section for the process 
e7 — > l + H depends on two free parameters m and m#. In our numerical estimation, 
we take the cross section o\ as a function of the lightest neutrino mass m for three values 
oirriH- Our numerical results are summarized in Fig. 2, Fig.3, and Fig.4, which correspond 
the subprocesses — > e + H , — >■ fi + H , and — > r + H , respectively. One can 
see from these figures that the value of the production cross section <j\ increases as m# 
increases, which is because the minimal value of the triplet VEV v& is proportional to the 
factor 1/rriH- In both of the NH and IH cases, the production cross section ai(fiH) for the 
subprocess —> ji + H is approximately equal to that for the subprocess cy — > t + H 
and their values are smaller than 16/6 in most of the parameter space of the HTM. 
This is due to h e/I w h eT , which reflects the neutrino mixing patterns. However, for the 
subprocess ey — > e + H , it is not this case. For the NH spectrum, its cross section 
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ai(eH) is larger or smaller than the cross section ai(fiH) ( or ax(rH)), which depends 
on the mass mo of the lightest neutrino. For the IH spectrum, the cross section a\(eH) 
is always larger than the cross section ai(^iH) ( or ui{tH)). For O.OOleV^ < itlq < leV 
and 300GeV < mu < QOOGeV, the values of the cross sections for production of H 
associated with a positron e + via collision at the ILC with yj~s = ITeV are in the 
ranges of 1.7/6 ~ 17.9/6 and 34.8/6 ~ 72.6/6 corresponding to the NH and IH cases, 
respectively. 
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Figure 5: The number of the e/i/i(a) and eee (b)events at the ILC. 



Ref. [9] has shown that, for < O.lMeV, the leptonic decays of iJ ±=t are the 
dominant decay modes, while the decay H ±ziz — > W^W^ 1 is negligible. In the context of 
the HTM, the leptonic decay channels if ±=l= — > l^l^ have been extensively studied in Refs. 
[11, 26] and the direct connections between their branching ratios and the triplet Yukawa 
couplings hij (i.e. the neutrino oscillation parameters) are discussed in these papers. For 
the leptonic decays of if ±=t , the possible decay modes are ee, fifi, tt, e/i, er and fir. In 
the case of 5 = and fi = <f2 = 0, /i/i and tt are the dominate decay modes, their 
branching ratios are roughly equal to each other for the NH spectrum and mi < O.leV. 
While for the IH mass spectrum and O.OOleV < m 3 < O.leV^, the doubly charged Higgs 
boson H mainly decays to ee and its branching ratio is 32% < Br(H — > ee) < 53%. 
Thus, for the final states generated by the subprocess ej — > e + H , we consider the three 
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leptons ~e\i\i and eee for the NH and IH cases, respectively. The number of the e/i/i and 
eee events at the ILC with y/s = ITeV and the integrated luminosity £i nt = 500/6 -1 
are shown in Fig. 5, where we have taken m in the range of O.OOleV ~ O.leV. From 
this figure, one can see that there will be several hundreds and up to ten thousands three 
lepton events to be generated per year in future ILC experiments. It is obvious that 
this is not the realistic number of the observed events. To take into account detector 
acceptance we should impose appropriate cuts. For example, the observed leptons must 
carry a minimum energy and respect a suitable rapidity cut, the angles of the observed 
leptons relative to the beam, 9 e , 6> M , must be restricted to some ranges [27]. Detailed 
analysis is needed, which is beyond the scope of this paper. 

It is well known that, at the ILC, the LFV processes can provide extremely clear 
signatures and are experimentally interesting. The leptonic decays of the doubly charged 
Higgs bosons give rise to same-sign lepton pairs, which can generate distinct experimental 
signals. The three lepton states e/i/i and eee generated by the subprocess ej — » e + H~ 
are almost free of the SM backgrounds at the ILC, which have been investigated in Ref. 
[27]. Thus, as long as the triplet VEV v& is enough small, the doubly charged Higgs 
boson H predicted by the HTM might be detected in future ILC experiments. 
4. Single production of the doubly charged Higgs boson H at the LHeC 

Recently, the high-energy ep collision has been considered at the LHC, which is called 
the LHeC [15]. For the LHeC, the energy E p of the incoming proton is given by the 
LHC beam, and the energy E e of the incoming electron is in the range of 50 — 200GeV, 
corresponding to the cm. energies of y/s = 2^E p E e 1.18 — 2.37TeV. Its anticipated 
integrated luminosity is at the order of 10 — lOO/fr -1 depending on the energy of the 
incoming electron and the design. Several studies have been performed to discuss the 
possibility of detecting the Higgs boson at the LHeC [28]. In this section, we consider 
single production of the doubly charged Higgs boson H via e^ collision at the LHeC . 

It is well known that the equivalent photon approximation (EPA) can be successfully 
to describe most of the processes involving photon exchange [29] . A significant fraction of 
pp collision at the LHC will involve quasi-real photon interactions occurring at energies 
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well beyond the electroweak energy scale. The quasi- real photons emitted by the incoming 
protons have low virtuality Q 2 and scattered with small angles from the beam pipe. The 
photon spectrum depending on virtuality Q 2 and energy E 1 can be described by the 
following relation [30, 31] 



with 



2 rn 2 p E 2 4m 2 p G 2 E + Q 2 G 2 M 



E P (E P -E 7 y * Am 2 p + Q 2 

G\ = % = (1 + %)\ F M = G 2 M , Q 2 = 0.71GeV 2 . (19) 

Where E p is the energy of the incoming proton which is related to the photon energy by 
E 1 = C,E P and m p is the mass of the proton. Thus, £ is the proton momentum fraction 
carried by the photon. The magnetic moment of the proton is /i 2 = 7.78, F E and F M 
are functions of the electric and magnetic form factors. Electromagnetic form factors are 
steeply falling as a function of Q 2 . Then, the Q 2 integrated photon flux can be written as 



4 max 



'W = jt"s^' (20) 

where Q 2 max ~ 2 — AGeV 2 . Since the contribution to the above integral formula is very 
small for Q 2 max > 2GeV 2 , in our numerical estimation, we will take Q 2 max ~ 2GeV 2 . 

At the LHeC, the effective production cross section a 2 (lH) for the subprocess cy —> 
l + H can be written as 



P^max 

a 2 (lH)= a(s)f(^M, (21) 

J (mii+mi Y s 



'(m/f +m;) 2 /s 

where s = AE e E^ = £s with E e = 140GeV and E p = 7TeV. To detect the intact proton 
or its loss energy, the forward proton detectors are needed. The ATLAS forward physics 
(AFP) collaboration will have forward detectors with an acceptance of 0.0015 < £ < 0.15 
[31]. The acceptance of the CMS-TOTEM forward detectors can reach 0.0015 < £ < 0.5 
[32]. In our numerical estimation, we will take £, max = 0.5 and £ m j n = (rriH + mi) 2 /s, which 
corresponds the cm. energy of the subprocess cy — > l + H , v 7 ? > m# + mi. 
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Our numerical results are summarized in Fig. 6 and Fig. 7, which plot the production 
cross sections for the subprocesses cy — > e + H and e^ — > [i + H at the LHeC as 
functions of the lightest neutrino mass mo for three values of the mass m# of the doubly 
charged Higgs boson. In these figures we have shown our results for the NH and IH 
cases. Our calculation shows that, at the LHeC , the effective production cross section for 
the subprocess —> fi + H is roughly equal to that for the subprocess — > r + H , 
which is similar with the conclusion at the ILC. Thus, we have not given the contours for 
the subprocess ey — > t + H . From these figures, one can see that the production cross 
section of the subprocess — > l + H at the LHeC is smaller than that for the ILC. 
For O.OOleK < m < leV and 300GeV < m H < 600GeV ', the values of the production 
cross section <j2{eH) and ^(/xi?) are in the ranges of 0.048/6 ~ 0.51/6 and 0.095/6 ~ 
0.43/6, respectively, for the NH spectrum, while are in the ranges of 0.97 fb ~ 2.03/6 
and 0.11/6 ~ 0.45/6, respectively, for the IH spectrum. If we assume the integrated 
luminosity Li nt = 100/6 -1 of the LHeC with -y/i = 1.96TeV, which corresponds E e = 
lAOGeV and E p = 7TeV, then there will be several and up to hundreds e + H events to 
be generated per year. 

The process ep — > e + H + X with H decaying to ee can give rise to number 
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of signal events with same-sign electron pair and an isolated positron, which is almost 
free of the SM backgrounds. Thus, it is also possible to detect the signatures of the 
doubly charged Higgs boson H via its production associated with a positron at the 
LHeC. Certainly, detailed confirmation of the observability of the signals generated by 
the HTM at the LHeC experiments would require Monte-Carlo simulations of the signals 
and backgrounds, which is beyond the scope of this paper. 
5. Conclusions and discussions 

Doubly charged Higgs bosons (i? ±=t ) appear in some popular new physics models 
beyond the SM, which can accommodate neutrino masses. This kind of new particles 
have distinct experimental signals through their decay to same-sign lepton pairs. Their 
observation in future high energy collider experiments would be a clear evidence of new 
physics. Thus, searching for H ±ziz is one of the main goals of current and future high 
energy collider experiments. 

The HTM is one of the attractive new physics models, in which a Higgs triplet is 
introduced and the neutrino masses can be obtained by the product of the triple VEV 
u& and the triplet Yukawa coupling hij. This model predicts the existence of the doubly 
charged Higgs bosons (if ±=t ). In this paper, we consider production of H associated 
with a positive lepton l + via cy collision at the ILC and LHeC. The production cross 
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sections are dependent on the neutrino oscillation parameters, the triplet VEV and 
the mass parameter m#. In our numerical estimations, we use the lower bound on v\ 
given by the LFV process /i — > ey, assume 5 = and <pi — <p 2 — 0, and take the 
experimental measurement values for the parameters s^, c^, and m^. In this case, the 
production cross sections depend only on two free parameters, the mass parameter m H 
and the lightest neutrino mass itlq. From our numerical results, we can obtain following 
conclusions. 

(1) The cross section for production of the doubly charged Higgs boson H associ- 
ated with a positive lepton l + at the ILC is generally larger than that at the LHeC. The 
production cross sections for the IH mass spectrum are larger than those for the NH 
mass spectrum at these two kinds of collider experiments. 

(2) Because of h epi ~ h eT , the production cross section of the subprocess — > fi + H~ 
is roughly equal to that of the subprocess cy — > t + H at both the ILC and LHeC 
experiments. For the IH mass spectrum, the cross section a(eH) is always larger than 
a(/iH), while for the NH mass spectrum, a(eH) is larger or smaller than <j({iH) which 
depends on the value of the lightest neutrino mass itlq. 

(3) For O.OOleV < m < leV and 300GeV < m H < 600GeV, the values of the cross 

section a(eH) at the ILC with = ITeV are in the ranges of 1.7/6 17.9/6 and 

34.8/6 72.6/6 corresponding with the NH and IH cases, respectively. While their 

values at the LHeC with y/s = 1.97TeV are in the ranges of 0.048/6 0.51/6 and 

0.97/6 2.03/6 for the NH and IH cases, respectively. 

(4) For the leptonic decays of H , the process ej — > l + H can produce three lepton 
final state, which is almost free of the SM backgrounds. In the context of the HTM, 
this process can only generate number of three lepton events at the LHeC , while, with 
reasonable values of the relevant free parameters, it can give rise to large number of three 
lepton events at the ILC. For instance, for O.OOleV < m < O.leV and 300GeV < 
rriH < 600GeV, the subprocess ej — > l + H with H decaying to ee can generate 
8.0 x 10 3 - 1.17 x 10 4 eee events at the ILC experiment with y/s = lOOOGeV and the 
integrated luminosity £ int = 500/6 -1 in the IH mass spectrum. Thus, the possible 
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signals of H predicted the HTM might be detected via this process in the future ILC 
experiments. 

In our numerical estimation, we have taken the minimal value of the triplet VEV 
z/a demanded by the current upper experimental bound for the LFV process \i — > ej. 
However, the constraint on the HTM from the process fi — > eee is generally stronger than 
that from the process \x — > e^y. For example, Ref. [12] has shown that, for tpi = if2 = 
and sin 2 2#i 3 = 0, the processes /i — > eee and /i — > cy demand v^mu > 400eV ■ GeV and 
VA m H > lOOeV • GeV, respectively. It is obvious that the minimal value of z/a demanded 
by /i — >• e7 is smaller than that for ji — > eee. Thus, if we use the constraints on the free 
parameters of the HTM coming from the process fi — > eee, the production cross sections 
and the number of the signal events will be decreased at least by one order of magnitude. 
However, as long as z/a < !0~ 5 GeV, there will be several eee events to be generated at 
the ILC. 

It has been shown [26] that the Majorana phases </?i and y2 2 have large effects on the 
branching ratios Br(H ±d: — > if if). Thus, </?i and y? 2 also have effects on the cross section 
for production of the doubly charged Higgs boson H associated with a positive lepton 
/ + . If we vary the values of <pi and <^ 2 , the numerical results for the cross sections a(eH), 
a{iiH), and g{tH) are also changed. However, our physical conclusions are not changed. 

Certainly, the doubly charged Higgs boson H ++ can also be singly produced at the 
ILC and the LHeC via the charge-conjugation processes of the corresponding processes 
for H . Similar to above calculation, we can give the values of the production cross 
sections for H ++ . Thus, the conclusions for H also apply to H ++ . 

If the doubly charged Higgs bosons H ±i: are indeed found at the LHC, various 
consistency checks will have to be performed. One of the important tasks is to exact 
calculate and precise measure the branching ratios Br(H ±i: — > if if) and see whether 
their values are consistent with a neutrino mass matrix from oscillation data. We expect 
that our work will be helpful to test the Higgs triplet mechanism for neutrino masses and 
further to distinguish different new physics models. 
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